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Abstract: Relatively little is known about fire regimes in grassland and cropland in Central Asia. In this 
study, eleven variables of fire regimes were measured from 2001 to 2019 by utilizing the burned area and 
active fire product, which was obtained and processed from the GEE (Google Earth Engine) platform, to 
describe the incidence, inter-annual variability, peak month and size of fire in four land cover types (forest, 
grassland, cropland and bare land). Then all variables were clustered to define clusters of fire regimes with 
unique fire attributes using the K-means algorithm. Results showed that Kazakhstan (KAZ) was the most 
affected by fire in Central Asia. Fire regimes in cropland in KAZ had the frequent, large and intense 
characters, which covered large burned areas and had a long duration. Fires in grassland mainly occurred 
in central KAZ and had the small scale and high-intensity characters with different quarterly frequencies. 
Fires in forest were mainly distributed in northern KAZ and eastern KAZ. Although fires in grassland 
underwent a shift from more to less frequent from 2001 to 2019 in Central Asia, vigilance is needed 
because most fires in grassland occur suddenly and cause harm to humans and livestock. 


Keywords: fire regime; burned area; active fire; K-means algorithm; Central Asia 


1 Introduction 


Fire regimes can be generalized to a particular combination of fire characteristics, including the 
frequency, intensity, variability, seasonality, size, extent, duration and pattern of fire (Krebs et al., 
2010). With global climate change, the frequent occurrence of extremely dry weather and the 
unreasonable use of fire by humans, and the instability of the factors influencing fire regimes is 
increasing (Flannigan et al., 2000; Flannigan et al., 2006; Moreira et al., 2011). Fire regimes are 
becoming the topics of worldwide concern because the occurrence of fire causes a series of social, 
environmental and economic problems (Pausas and Fernandez-Mufioz, 2012; Chuvieco et al., 
2014). Fire also poses serious threats in Central Asia, especially in Kazakhstan (KAZ). To better 
understand and manage fire regimes, there is an increasing interest in monitoring fire distribution 
and assessing potential drivers (Calvifio-Cancela et al., 2017). However, to date, large-scale 
schemes for fire regimes are limited, and few studies have considered the impact of fire on 
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grassland in Central Asia. 

Data for mapping fire regimes can be obtained from historical records (interviews, diaries and 
literature searches), palaeoecological evidence (charcoal in lakes, soil sediments and tree scars), 
and remote sensing data of fire patterns (Morgan et al., 2001). Considering the integrity, 
availability, spatial convergence and relative cost of data records, satellite data are often 
employed for reconstructing maps of fire regimes (Downing et al., 2017). In recent years, as the 
algorithms for burned area of fires, including MCD45A1, MCD64A1 and FireCCI5.1, have 
become more powerful and been widely calibrated (Justice et al., 2002), products have been used 
to map fire regimes in real time and with spatial extents (Morgan et al., 2001; Levin and 
Heimowitz, 2012; Chen et al., 2017), and extensive research on biomass burning at a spatial 
resolution of 500 m and approximately 250 m has been conducted (Archibald et al., 2009; 
Petrenko et al., 2012; Vadrevu et al., 2012; Grégoire et al., 2013; Lehmann et al., 2014; Oliveras 
et al., 2014; Yang et al., 2014; Hantson et al., 2015; Chen et al., 2016; Andela et al., 2017). Active 
fire represents the approximate central site of a fire, but the shape and size of a large fire cannot 
be accurately portrayed (Roy et al., 2013). Therefore, burned area products play a key role in 
mapping medium and large fire regimes. To address the limitations of these products, researchers 
applied active fire and burned area data to map fire regimes in numerous studies (Devineau et al., 
2010; Chen et al., 2017). The Landsat satellite has a high spatial resolution, which allows fire 
regimes to be monitored using spectral indices, including the normalized burn ratio (NBR) index, 
burned area index (BAI), mid-infrared burn index and global environment monitoring index 
(GEMI). And the effects of a single fire analyzed by carbonization after burning and the 
disappearance of plants or by analyzing ecological recovery across a region was estimated 
(Chuvieco et al., 2002; Bastarrika et al., 2011; Schepers et al., 2014; Chen et al., 2017; Downing 
et al., 2017; Musyimi et al., 2017; Argibay et al., 2020). However, Landsat images are typically 
used to analyze the differences between pre- and post-fire images in a specific research region. 
Due to the contamination of Landsat cloud images and a limited time scale that cannot distinguish 
between new and old fire scars, MODIS (Moderate Resolution Imaging Spectroradiometer) is 
more popular method for characterizing the temporal and spatial distribution of fire regimes in 
very broad geographical regions (CHUVIECO et al., 2008). 

As fire is rarely studied in Central Asia, the results of this study will aid in filling this research 
gap. In this study, we analyzed all the available time series images (FireCCI5.1, Fire Information 
for Resource Management System (active fire data) of Central Asia in Google Earth Engine 
(GEE)) (1) to develop an easy and robust method to map fire regimes; (2) to utilize this method to 
reconstruct and characterize maps of fire regimes in Central Asia; (3) to analyze the monthly 
number of active fires in different land cover types (grassland and cropland) from 2001 to 2019 in 
Central Asia; and (4) to discuss the impact of fire on grassland productivity, the effects of drought 
on fires in grassland, and seasonal differences of fires in cropland in Central Asia. 


2 Materials and methods 


2.1 Study area 


Central Asia (35°07'-55°26’N, 46°29'-89°19’E) is located in the Eurasian mainland and has a 
semi-arid and arid climate (Qi et al., 2012). The study area includes five countries (KAZ, 
Kyrgyzstan (KGZ), Tajikistan (TJK), Turkmenistan (TKM) and Uzbekistan (UZB)), which are 
collectively referred to as Central Asia in this study. As shown by the MODIS land cover product 
in 2018 (Fig. 1), approximately 75% of the land in Central Asia is covered by grassland, but the 
coverage density is uneven based on the mean value of leaf area index (LAI). Forest is mainly 
located in eastern KAZ, northern KAZ and Aqmola in KAZ. According to an analysis of the 
annual forest area from 2001 to 2018 based on the MODIS MCD12Q1 V6 product, it is surprising 
to find that the area of forest in KAZ has undergone a process of large-scale deforestation and 
gradual recovery. Agriculture can be roughly divided into two types in Central Asia. Rainfed 
cropland is mainly distributed in the northern KAZ, including North KAZ, Aqmola and Qostanay 
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in KAZ, which mainly grows wheat and barley, and this region is also one of the most important 
wheat exporters in the world. Irrigated cropland is mainly located in central and southern Central 
Asia, (Jiang et al., 2017). 

As indicated by the long-term climate record, the climate of Central Asia is generally arid, 
especially in the southern UZB and northern TKM, where the average annual precipitation is less 
than 100 mm, but there are relatively humid regions, such as the northern KAZ (approximately 
400 mm) and eastern mountainous regions of Central Asia (600-800 mm) (Klein et al., 2012). In 
recent decades, Central Asia has experienced severe and growing water deficits with increasing 
temperature and variable precipitation patterns, causing frequent drought events (Guo et al., 2015; 
Guo et al., 2017). Simultaneously, unmanaged fires in cropland put great pressure on firefighting 
in Central Asia. 
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Fig. 1 Main types of land use in Central Asia in 2018 (a) and fire burn scars based on Landsat 8 images in 
central KAZ (b—-g, marked as A-F in Figure 1a). A, northeastern Aktobe; B, northeastern Aktobe; C, southern 
Qostanay; D, southern Qaraghandy; E, southern Qaraghandy; F, East Kazakhstan. KAZ, Kazakhstan; UZB, 
Uzbekistan; TKM, Turkmenistan; KGZ, Kyrgyzstan; TJK, Tajikistan. Abbreviations are the same as in the 
following figures. 


2.2 Data sources 


Three MODIS products were used in this study: (1) MODIS active fire data originating from the 
Fire Information for Resource Management System; (2) burned area FireCCI5.1 product; and (3) 
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MODIS land cover type product collection 6 (MCD12Q1 V6). 

The Fire Information for Resource Management System (FIRMS) dataset contains the active fire 
product processed by Land Atmosphere Near Real-time Capability for EOS (LANCE) that 
integrates MOD14/MYD14 fire and thermal anomaly products in raster form, and the active fire 
location represents the centroid of a 1-km pixel (Davies et al., 2009). 

The FireCCI5.1 products, which originated as part of the European Space Agency (ESA) Climate 
Change Initiative (CCI) program, are new monthly global burned area products at a resolution of 
approximately 250 m (Lizundia-Loiola et al., 2020). To analyze the characteristics of fire regimes 
including the duration, size, area, variability and seasonality, we processed burned area (Fire CCI 
5.1) and active fire (FIRMS) products from 2001 to 2019. 

The MCD12Q1 collection from six land cover products (MCD12Q1 V6), which is generated 
using supervised classifications of MODIS Terra and Aqua reflectance data, provides annually 
updated global land cover data from 2001 to 2019 at a spatial resolution of 500 m (Zeng et al., 
2020). The MCD12Q1 V6 improved upon those previous versions and contained five global land 
classification systems. The International Geosphere-Biosphere Program (IGBP) classification 
scheme is selected in this study. The seventeen land cover classes from the IGBP classification 
scheme were merged into four types of land cover with vegetation: forest, grassland, cropland and 
bare land. Non-vegetation land cover types containing permanent wetland, urban areas, built-up 
areas, permanent snow, ice and water body were precluded from the study. These variables were 
obtained from GEE and resampled to a resolution of 250 m. 


2.3 Methods 


Eleven variables were measured to characterize the fire regimes and their attributes that roughly 
reflect the time, space, size and peak month of fires in Central Asia (Table 1). These attributes 
consider the annual average characteristics of fires, inter-annual variability, size distribution, fire 
duration, peak month and effects of vegetation. Variables 1, 3 and 5 were calculated using burned 
area data, while variables 2, 4, 6, 7, 8, 9, 10 and 11 were calculated using active fire data. 


Table 1 List of fire variables from 2001 to 2019 


No. Type Abbreviation Variable Unit 
1 ABAY Average annual burned area hm?/a 
Fire incidence 
AFOF Average annual active fire frequency counts/a 
3 CVBA Inter-annual coefficient of variance in annual burned area 
Inter-annual variability oa 2 : ceed 
Inter-annual coefficient of variance in annual active fire 
4 CVFP 
frequency 
5 Fire size GI Gini index 
6 Active fire FPSD Seasonal duration of fires d 
7 seasonality FPDT Peak month of fires 
8 PFAF Percentage of forest active fires 
9 PGAF Percentage of grassland active fires 
Vegetation type % 
10 PCAF Percentage of cropland active fires 2 
11 PBAF Percentage of bare land active fires 


The Gini index (GI) concerns fire size and characterizes the fire size distribution tendency. In 
other words, GI is applied to characterize the extent to which the burned area is concentrated on a 
small number of larger fires or is recognized as a large number of smaller fires rather than a single 
large fire (Chen et al., 2017). A higher GI value indicates that the burned area is concentrated in 
some large fire regions. The duration of active fire involves monitoring the duration of the period 
from the occurrence to the duration of the fire using the confidence band of FIRMS dataset from 
10% to 90% (Chen et al., 2017). We calculated the peak month of fires based on the month in which 
the total number of active fires in each month from 2001 to 2019 was the largest in a specific 
location (Chen et al., 2017). The percentage of each of the four types of vegetation affected by fires 
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is based on the ratio of the active fire counts in each vegetation type to the total active fire counts in 
the grid cell (Chen et al., 2017). 

The eleven calculated indicator maps were integrated into a multiband image in GEE. The tool 
for the unsupervised classification K-means algorithm in GEE will be applied to classify images 
using 5000 sample points. Computing the Euclidean distance between the sample points allows the 
smaller similarity and discrepancy to be gathered into one class, and multiple clusters are eventually 
formed so that the sample within the same cluster has high similarities and large differences among 
different clusters. The advantages of this classification method are its simple principle, easy 
implementation, fast convergence speed and obvious differences among classifications (Steinley, 
2006), which is very effective for processing high-dimensional images. Moreover, this method 
effectively conserves the original cluster structure (Steinley, 2006). Therefore, K-means algorithm is 
applied in this study. 

In terms of the given k clusters, the K-means algorithm classifies the given data. However, 
choosing the right number of clusters is very important for classifying results. A high k value does 
not necessarily signify a better quality of information (Benmouiza and Cheknane, 2013). On the 
contrary, a small number of clusters can cause ambiguous results that do not convey the true 
meaning. Based on the 5000 sample points obtained from the image, we determined the optimal 
number of clusters through analysis to obtain the most robust K-means results. The clusGap 
function in the R language cluster function package uses Gap statistics to determine the number of k 
(Furihata et al., 2019) and uses bootstrap sampling of the data to compare internal differences (500 
samples are selected in this article). This method is considered more effective for determining the 
value of k than other methods (Tibshirani et al., 2001). The results of the evaluation are mainly 
based on silhouette analysis and the number of cluster statistics according to Gap statistics. 


3 Results 


3.1 Characteristics of fire regime variables 


Figures 2 and 3 showed the annual fire incidence in Central Asia and its respective inter-annual 
coefficient of variance from the burned area products and active fire products. In Figure 2, the 
ABAY reveals that the largest burned area occurred in KAZ and the Nukus Irrigation District in 
UZB. In addition, ABAY also occurred in the Ahal and Mary regions of TKM. KAZ was the 
country with the worst fires in Central Asia, especially in the eight states in the north-central 
region of KAZ (West KAZ, Aqtobe, Qostanay, North KAZ, Aqmola, Pavlodar, Qaraghandy and 
East KAZ) and the Syr Darya and Chuhe River basins. Similarly, Figure 3 shows AFOF in each 
year, confirming the frequent fires in the abovementioned areas of Central Asia. In TJK and KGZ, 
the area and frequency of fires were low. It is possible that the fire burn algorithm of FireCCI5.1 
has a high threshold in the mountainous area and could not detect an active fire or burned area 
there. 

As expected, in the regions with high fire occurrences i.e., ABAY and AFOF, the observed 
inter-annual coefficients of variance are generally lower, and in regions with low fire occurrences, 
the stability of the coefficient of inter-annual variability is weak, especially in the inter-annual 
coefficient of variance in active fires, i.e., CVFP. The inter-annual coefficients obtained from 
different datasets in the central region of KAZ were found to show a strong variability. The 
coefficient of variation around the northern region of KAZ is small, and the same phenomenon is 
also found in the Chuhe River Basin in KAZ and Nukus Irrigation Districts in UZB and in the 
Ahar and Mary regions of TKM. 

In Figure 4, the FPDT occurred in spring in northern KAZ (Qostanay, North KAZ and Aqmola), 
the Syr Darya and Chuhe River basins in KAZ, and the Nukus Irrigation District in UZB. In 
summer, FPDT was mainly distributed in the central region of KAZ, including West KAZ, Aqtobe, 
Qostanay, Qaraghandy and East KAZ. Moreover, Tashkent city and the Fergana Valley in KGZ 
had a high-density region of fire. In TKM, FPDT was unevenly distributed and mainly occurred 
in winter and summer seasons. Overall, depending on the month, the locations of high-density 
fires varied. 
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Fig. 2 Average annual burned area (ABAY, a) and its inter-annual coefficient of variance (CVBA, b) in Central 
Asia based on FireCCI5.1 burned area data 
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Fig. 3 Average annual active fire frequency (AFOF, a) and its inter-annual coefficient of variance (CVFP, b) in 
Central Asia based on MODIS active fire data 
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Fig. 4 Peak month of fires (FPDT) in different months in Central Asia based on MODIS active fire data 


FPSD result is presented in Figure 5. The regions with the longest durations in KAZ include 
North KAZ, Qostanay, Aqmola, East KAZ, Chuhe River Basin, Syr River Basin and the region 
bordering KAZ and KGZ. The regions with the longest durations include the Nukus Irrigation 
District and near Tashkent city in UZB. The longest durations are in the states of Ahal and Mary 
in TKM. Active fires in central KAZ have shorter durations. In Figure 5, GI shows that in KAZ 
region, especially in the northern part of Qostanay, North KAZ, Aqmola, Chuhe River Basin, Syr 
River Basin and the border region of KAZ and KGZ, GI value is high. In addition, GI is high in 
the Nukus Irrigation District of UZB. In central KAZ, GI is low, indicating that there are many 
small fires and that the burned area is more evenly distributed or that there are few fires in the 
region. In addition, GI value in parts of KGZ and TJK was high this time, indicating that fires in 
the region were concentrated and that there were a few large fires. 
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Fig. 5 Seasonal duration of fires (FPSD, a) and Gini index (GI, b) in Central Asia based on MODIS burned area 
data and active fire data 


The vegetation type affected by fire was obtained from the annual MODIS land classification 
IGBP product. The forest distributed in Central Asia is relatively small and less affected by fires, 
and major fires were found in the states of East KAZ and North KAZ (Fig. 6). The vegetation 
type most affected by fires in Central Asia was grassland, especially in KAZ, as well as in the 
Nukus Irrigation District and the Ahar and Mary regions in TKM. Farmland burning and fires 
occurred at about the same time in the region. And it was found that the frequency of farmland 
burning was generally high. The occurrence of farmland burning is extremely frequent in Central 
Asia. Moreover, fires in cropland differed in the time of burning straw due to the ripening 
variation (Fig. 6). 
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Fig. 6 Vegetation types affected by fires in Central Asia based on MODIS active fire data. PFAF, percentage of 
forest active fires; PGAF, percentage of grassland active fires; PCAF, percentage of cropland active fires; PBAF, 
percentage of bare land active fires. 


3.2 Clustering of fire regimes 


Figure 7 shows that fires in Central Asia present an unbalanced spatial distribution, and many 
indicators have outliers. The figure shows that three variables (ABAY, AFOF and GI) reached the 
highest or the lowest values in few locations, indicating that fire burned areas, frequency and fire 
size in few regions were much larger or smaller than those in other regions in Central Asia. 

The K-means algorithm is classified by selecting an appropriate number of clusters. Gap 
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statistics are calculated by the clusGap function to obtain the best classification number. A higher 
silhouette value is considered to be one of the prerequisites for successful clustering (Benmouiza 
and Cheknane, 2013). A 0.60 silhouette value is a good result for clusters (Lleti et al., 2004). 
Surprisingly, the silhouette value obtained in this experiment is 0.88, which meets the 
requirements of experimental analysis. After silhouette analysis, k=5 is the best classification. 
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Fig. 7 Box plot distribution of fire variable values. The black center mark indicates the median, while the 
bottom and top edges of the green box and line indicate the 25'" and 75" percentiles, respectively. Outliers (if any) 
are indicated by red plus signs. 


In terms of the median tendency of eleven variables, we divided the characteristics of the five 
fire clusters into three to four quantitative levels, and the results are summarized in Table 2. 
Finally, a map of fire regimes in Central Asia is given in Figure 8. It is concluded that 43% of the 
land in Central Asia is threatened by fires. Table 2 lists five types of fire regimes that are different 
from the above clustering process. The five types of fire regimes are described as follows: 

Fire regime 1: Mainly distributed in bare land, these areas are less threatened by fires and have 
low fire frequency, which may be affected by the low-density fuel bed, accounting for 57.0% of 
the overall fire-affected area. 

Fire regime 2: This regime represents fires in cropland, including medium burned areas, 
frequent active fires, high inter-annual variabilities in burned areas, large intense fires and low 
fire duration, occupying 5.8% of the overall fire-affected area in the southwest of Central Asia. 
The peak months of fires are mainly in winter and summer, which may be caused by crop 
ripening and farming habits, which are discussed later. 

Fire regime 3: This regime is characterized by a large burned area, a medium frequency and 
inter-annual variability of active fire, and by more frequent long duration of fires in grassland in 
July, and it mainly occurs in central KAZ, with 17.9% of the overall fire-affected area. 

Fire regime 4: This regime is characterized by an irregular burned area, low frequency small 
fires, and more fires in grassland in July, with a high inter-annual variability, accounting for 
15.8% of the overall fire-affected area. 

Fire regime 5: This regime is mainly distributed in northern KAZ, the Syr River Basin and 
Nukus Irrigation District, which belongs to the agricultural area in Central Asia. This regime is 
characterized by large burned areas and a high frequency of fires that are large and intense, with 
long duration, and it accounts for approximately 4.0% of the overall fire-affected area. The peak 
months of fire mainly occur in spring. 
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Table 2 Fire regime variables and classification of fire regimes 


Fire regime 


Index 
1 2 3 4 5 
ABAY (hm?/a) Low (0.10) Medium (0.30) High (0.50) Medium (0.30) High (0.34) 
AFOF (counts/a) Low (0.10) High (0.20) Medium (0.15) Low (0.10) High (0.30) 
CVBA Low (2.00) High (4.20) Low (2.80) High (4.00) Medium (3.00) 
CVFP Low (1.70) Medium (3.00) Medium (2.80) High (4.20) Low (1.90) 
GI Low (0.71) High (0.85) Medium (0.82) Low (0.65) High (0.88) 
FPSD (d) 0.00 Low (0.45) High (3.00) Medium (2.00) High (5.00) 
FPDT - 1.00 7.00 7.00 4.00 
Vegetation type Bare land Cropland Grassland Grassland Cropland 


Note: -, no value; high-, medium- and low-grade medians are shown in brackets. 
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Fig. 8 Distribution of fire regime in Central Asia 
3.3 Time series of active fires in grassland and cropland 


As fires in grassland occur very suddenly and are destructive, time series analysis is expected to 
provide valuable information for fire prevention and control. The number of monthly active fires 
from 2001 to 2019 in Central Asia was used to analyze seasonal changes of fires in grassland and 
cropland. Figure 9 shows that the fire regime conditions across the countries exhibit different 
seasonal changes over time. Before 2010, prairie fires in KAZ were mainly concentrated in April 
to October, after which the number of active prairie fires suddenly decreased. Before 2008, KGZ 
had many active fires in April, August, September and October, after which the number of active 
fires gradually decentralized. The year 2007 was an important turning point for active fire 
distribution. Before 2007, the active fires in TJK mostly occurred in March, and from June to 
November, the fires in grassland in UZB were concentrated from March to October, while for 
TKM, the active fires were mainly distributed in the two periods of February to March and June 
to October. After 2007, the active fires were concentrated in March, August and September for 
TJK, March to April for UZB, and February to March for TKM. In general, the active fires in 
grassland significantly decreased in recent years in Central Asia. 

Figure 10 shows that in Central Asia, all the countries used seasonal straw burning to improve 
soil quality. The months of straw burning in KAZ are concentrated between April and May. For 
KGZ, the months of straw burning are March to April and August to October, and the straw 
burning period in TKM is similar to that in KGZ, concentrated in the periods from February to 
March and June to July. In UZB, this practice takes place in July. Different from the other 
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countries, the months of straw burning in TJK is relatively scattered. However, straw burning not 
only pollutes the air but also causes large-scale fires in grassland. 
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Fig. 9 Monthly active fires in grassland in Central Asia 


4 Discussion 


4.1 Impact of fire on grassland productivity 


Fire cannot only quickly reduce or remove the biomass on the surface, but it does increase the soil 
temperature, resulting in a decrease in the soil water content and affecting the water use efficiency 
of vegetation (Zavala et al., 2014). Some scholars have stated that fire promotes the renewal of 
grassland ecosystems and even produce dependence (Pereira et al., 2013). In the central region of 
KAZ, grasslands are frequently burned by fire. To explore the impact of fire on the productivity 
of grassland in the region, we selected the annual MODIS net primary productivity (NPP) product 
of vegetation (MOD17A3HGF V6) to study how the NPP of vegetation varies after the grassland 
has been burned by fires. 
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Fig. 10 Monthly active fires in cropland in Central Asia 


All six regions including northeastern Aktobe (A and B), southern Qostanay (C), southern 
Qaraghandy (D), southern Qaraghandy (E) and East KAZ (F) were burned from August to 
September 2016 (Fig. 11). Of these regions, A, B and C belong to Fire Regime 3 and are 
characterized by a large burned area, high-intensity and long-duration fires. D, E and F belong to 
Fire Regime 4 and are characterized by an irregular burned area and small low-intensity fires. A, 
B and C were burned in 2016, and the results show that the grassland NPP began to decline in 
2017, which reflects the poor recovery ability of vegetation. However, for Fire Regime 4, the NPP 
of grassland began to decrease in 2016 after burning and then gradually increased, which 
indicates that moderate-scale fires are conducive to regeneration and restoration in grassland. Of 
course, the influencing factors of grassland productivity have a considerable relationship with 
precipitation, temperature, grazing and human activities (Stavi, 2019). Nevertheless, in the long 
run, Fire Regime 4 leads to the better vegetation recovery ability than Fire Regime 3 after a fire. 


4.2 Seasonal differences of fires in grasslands across countries 
In Central Asia, active agricultural fires are closely related to straw burning. Due to the large 
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Fig. 11 Net primary productivity (NPP) of vegetation in grassland from 2016 to 2019 based on different fire 
regimes. A, northeastern Aktobe (25 August, 2016); B, northeastern Aktobe (15 August, 2016); C, southern 
Qostanay (4 September, 2016); D, southern Qaraghandy (21 August, 2016); E, southern Qaraghandy (11 August, 
2016); F, East Kazakhstan (1 September, 2016). 


volume of straw and high transportation costs, burning is the least expensive method of removing 
straw to prepare land after harvesting crops (Long et al, 2016). The burning time of straw is 
closely related to the local planting and harvesting times. 

Normalized difference vegetation index (NDVI) in agricultural areas reflects a phenological 
cycle of crops from growth to ripening (Zhang et al., 2019). In other words, there is a crest in the 
time series, which means that the crops in the region ripen once a year. If there are two crests, it 
means that the region has two ripening events once a year. If there are three peaks, it means that 
there are three ripening events once a year (Zhang et al., 2019). To explore the reasons for the 
seasonal differences of fires in cropland in Central Asia, we obtained long-term NDVI products 
(MOD13Q1 V6) from GEE, and drawn NDVI change trends on the platform in major agricultural 
regions. 

In Figure 12, there is a peak in KAZ and KGZ, indicating that the main agricultural regions of 
the countries ripen once a year, which explains why the time of straw burning of both countries 
occurred in spring and autumn. In UZB, there is a peak in the Nukus Irrigation District, and there 
are two peaks in southern Tashkent and the Fergana Valley, which shows that the land will be 
replanted in summer when straw burning occurs. In UZB, the agricultural areas where can be 
planted twice a year are larger than that of Nukus Irrigation District, where can be planted once a 
year, and more fires occur in summer. The agricultural areas of TKM and TJK planted crops twice 
a year and fires in cropland occur in summer twice too. 


4.3 Effects of drought on fires in grassland 


Drought is an important factor that causes frequent natural fires. In drought-prone regions, such as 
Central Asia (Guo et al., 2018a, b, 2019), due to long-term high temperatures and scarce 
precipitation, the vegetation dries quickly, which makes these areas more prone to fires. To 
explore whether or not drought severity is more likely to lead to a fire, we further discussed the 
relationship between fire regime and drought conditions in grassland in KAZ at the annual and 
monthly scales. The Palmer Drought Severity Index (PDSI) (Palmer, 1965) was used to judge the 
drought condition because it can successfully quantify the severity of droughts across different 
climates (Wang et al., 2015). 

Figure 13a show the PDSI (Palmer Drought Severity Index) drought index change from 2001 to 
2018, which indicates that the ratios of near normal years, wet years and drought years are 0.39, 
0.22 and 0.39, respectively. The values suggest that Central Asia was dominated by drought from 
2001 to 2018. Does a higher degree of drought mean a higher risk of fire? We further analyzed the 
correlation between PDSI of annual fire accumulation frequency and cumulative area (Fig. 13b), 
and found that the annual cumulative burned area and fire frequency in grassland do not correlate 
well with drought. Therefore, at the annual scale, the drought condition is not the main factor 
affecting the cumulative burned area and frequency of fires in grassland in Central Asia. 
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Fig. 12 Time series of normalized difference vegetation index (NDVI) in the agricultural areas of Central Asia 
from 2015 to 2019 based on MOD13Q1 V6 products. KAZ1, northern Kazakhstan; KAZ2, Syr Darya Valley; 
KGZ1, northern Kyrgyzstan; KGZ2, Uzgen City; UZB1, Nukus Irrigation District; UZB2, southern Tashkent City; 
UZB3, Fergana Valley. 
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Fig. 13 Relationship between PSDI (Palmer Drought Severity Index) and annual burned area in grassland. (a), 
PDSI change from 2001 to 2018; (b), correlations of PDSI with frequency and cumulative area of fires in 
grassland in Kazakhstan. 


The relationship between burned area in grassland and PSDI indicates that near normal and 
humid conditions occurred in most months, while the distribution of dry months gradually 
changed from random to concentrated. From Figure 14a, we can see that drought was mainly 
concentrated from April to October, especially in June and September, which is consistent with 
the high-frequency active fires in grassland (Fig. 14b). In winter, regardless of drought or wet 
conditions, the burned area in grassland was very small. In summer, persistent droughts often 
result in large-scale fires in grassland. According to statistics, when PDSI reveals a moderate or 
more severe drought event in summer, the probability that the burned area in grassland will 
exceed 5000 km? is 73.3%, which may be related to the highest vegetation biomass at the time. 
Moreover, it was found that the large monthly burned area in grassland (1x10*-3x10* km?) is 
negatively correlated with PDSI (r= —0.41, P=0.017) and strongly positively correlated with PDSI 
(r=0.82, P=0.020) after more than 35,000 km? monthly burned area in grassland. It can be seen 
that a certain scale of burned grassland is closely related to the drought at the monthly scale. 
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Fig. 14 Relationship between PDSI (Palmer Drought Severity Index) and monthly burned area in grassland. (a), 
monthly PDSI from 2001 to 2018; (b) monthly burned area in grassland from 2001 to 2018 in Kazakhstan. 


5 Conclusions 


The overall distribution of fires in Central Asia is uneven, and KAZ is the country most 
threatened by fires. Grassland in KAZ is the most affected vegetation type by fire. The fires in 
grassland in five countries have gradually decreased in recent years, which show that the fire 
control by government has achieved remarkable results. The fires in cropland are mainly caused 
by straw burning, and the burning time of straws in different countries is inconsistent. However, 
due to the suddenly occurrence of fires in grassland and huge effects on humans and livestock, 
vigilance is needed. There are still some shortcomings in our study, especially the fact that 
ignition source as an important component of a fire regime was not involved in the work, because 
of the lack of ground fire survey data. Therefore, it is difficult to distinguish between natural fire 
sources and man-made fire sources only by remote sensing data. As a research topic to be 
explored in the next stage of work, a fire risk analysis will be carried out under the interactive 
influence of human activities and natural factors. Moreover, we will use a variety of indicators 
and algorithms to test the correlation between fire variables and ecological and human factors. 
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